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The fundamental solutions of diffusion equation for the local-equilibrium and nonlocal models are considered as 
the limiting cases of the solution of a problem related to consideration of the Brownian particles random walks. 
The differences between fundamental solutions were studied. It was shown that on the period of observation 
time and distances exceeding the time and space associated with one step of random walk the fundamental 
solutions of diffusion and telegraph equations are very close to each other. In particular, these fundamental 
solutions are very close to the values of the probability density distribution of diffusive particles that is 
obtained from the accurate solution of the random walk problem. The difference between the probability 
density values for the models considered is decreased with time as t -3 / 2 , the relative difference as t _1 , and 
velocity of propagation of these fixed distribution densities is proportional to t 1 / 2 . The new modified non-local 
diffusion equation is suggested. It contains only microparameters of the random walk problem. 
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I. INTRODUCTION 

It is known that one of the properties of linear diffusion 
equation solution 


dC(x, t) 

dt 


DAC(x,t ) 


(1) 


is the infinite velocity of perturbations propagation. It is 
seen, for example, from the fundamental solution of this 
equation 


Pc(x,t ) 


m 

V'lnDt 




( 2 ) 


connecting the current particle function and the concen¬ 
tration gradient. By analogy, for conductive heat trans¬ 
fer one can write the Fourier law connecting the heat 
flux with temperature gradient and for liquids and gases 
involved in the filtration phenomenon in porous medium 
the Darcy’s law connecting the stream of liquid with pres¬ 
sure gradient. 

One of the approaches related to resolution of the prob¬ 
lem connected to corrections for the infinite velocity of 
perturbations propagation lies in the usage of linear non¬ 
local models, for example, suggested by the extended of 
nonequilibrium thermodynamics^. In this case in ex¬ 
pression (3) some additional relaxation terms are intro¬ 
duced. In the simplest case, the Fick’s law in the nonlocal 
interpretation is written as 


Here pc{x,t) is the fundamental solution presented in 
the Gaussian form, C is concentration, D is the diffusion 
coefficient, 0 defines the unit step Heaviside function, x is 
associated with current coordinate, t is the current time. 
This solution contradicts to intuitive speculations that 
diffusive particles and any information about their move¬ 
ment are propagated with finite velocities. We should 
mark that equation © is a consequence of supposition 
that principles associated with locality and local ther¬ 
modynamic equilibrium are valid. In the frame of this 
supposition one can use the diffusion Fick’s law in stan¬ 
dard form 


J(x,t) + r d = —DS7C(x,t) (4) 

and equation © is reduced to the telegraph equation of 
the type 


d(x,t) d 2 {x,t) 
dt +T dt' 2 


D A(x, t). 


( 5 ) 


The fundamental solution of the last one for ID case is 
expressed as 


J(x,t) = -DVC(x,t) 


( 3 ) 


Pte(x, t ) 


. ^ exp(——) 0 
y/l D 7 2 t’ 



h 



9 Electronic mail: Marat.Ovchinnikov@kpfu.ru 


1 

27 


( 6 ) 












2 


where Iq is the modified Bessel’s function of the first kind 
having zeroth order. The paradox of the infinite veloc¬ 
ity is eliminated and the limiting velocity describing the 
perturbations propagation in a nonequ ilibriu m medium 
becomes finite and is equaled to V = y/D/r. 

In the given paper the proper selection of the fun¬ 
damental solutions plays a key role. Their knowledge 
determines the specific solutions of linear equations of 
mathematical physics as well, because the correspond¬ 
ing Green’s functions can be expressed in the form of 
linear combinations of fundamental solutions including 
their derivatives and integrals. 

Equations of type CD for many real situations of 
macroscopic physics are good approximation. But to¬ 
gether with this statement the consideration of the non¬ 
local effects is an actual problem for relatively small inter¬ 
vals of observation times and distances. It is interesting 
to consider the evaluation of the limits of applicability of 
the local thermodynamic equilibrium hypothesis and the 
finding of approximate approaches related to determina¬ 
tion of the current of the value studied and gradient of 
potential of the corresponding scalar field. Taking into 
account the fact that approach ( 4 ) is also approximate, 
it would be interesting to evaluate errors which arise in 
the usage of similar expressions in comparison with solu¬ 
tion of the ’’exact” transfer equation which is remained 
unknown. 

The certain complexity in experimental testing in lab¬ 
oratory conditions the adequacy of equations © and 
( 4 ) lies in nonlinear dependencies of the key parameters 
involved. It creates some difficulties in interpretation of 
the results obtained. So, in the given paper the dynam¬ 
ics of hypothetical linear systems modelling the transfer 
phenomenon is considered. 

It is necessary to remind here about the inertia heat 
phenomenon when nonlinearity leads also to the finite ve¬ 
locities of perturbations^ but similar and relatively exotic 
models are out of the scope of this research. 

If one come back to micro-level the diffusion phe¬ 
nomenon, conductive heat transfer and filtration of liq¬ 
uids and gases in porous media and similar phenomena 
can be described with the help of random behavior of mi¬ 
croparticles. For diffusion and heat transfer phenomenon 
one can introduce concepts as the characteristic length 
Ax and time At of free particles path or quasiparticles of 
phonon type; for filtration we consider the characteristic 
distances and times associated with the hopping process 
of liquid droplets from one pore to another one. 

Hence, from microscopic point of view of the phe¬ 
nomena considered it is important to establish the cor¬ 
relations between phenomenological parameters of the 
corresponding theories similar to the extended irre¬ 
versible thermodynamics with characteristic relaxation 
times and microparameters of the transfer process stud¬ 
ied. 

In random walk model the hopping length Ax and 
hopping time At characterize the dynamics of the sys¬ 
tem considered on microlevel. For various media and 


processes these microparameters are varied in wide lim¬ 
its. The length of free path for phonons in crystals is 
changed in the range 10“ 7 — 10“ 10 m that corresponds 
to the free path time located in the range 10” 10 — 10“ 13 
s for frequencies of the phonon spectrum from the in¬ 
terval 10 10 — 10 12 Hz. We should note that in paper— 
with the usage of the molecular dynamics methods it 
was shown that characteristic relaxation times in nonlo¬ 
cal heat transfer model written in the form of telegraph 
equation ( 5 ) for solid argon crystal are very small ( 10 -11 
s) that comparable with times of the phonons free path in 
a crystal. The characteristic microscopic hopping length 
in diffusion phenomenon in the most cases is located in 
the range 10“ 3 — 10“ 9 m, for filtration 10“' — 10“ 3 m 
and the corresponding characteristic times are located in 
the intervals 10“ 10 —10° s and 10“ 3 —10 _1 s, accordingly. 


II. THE MODEL 

The random walk model at present time is widely used 
for description of many processes; in particular, the frac¬ 
tal dynamic models are developed in-. Supposing that 
on validity not only nice equations can pretend but also 
the nice algorithms, as well, we will be based on suppo¬ 
sition that true description of diffusion phenomenon is 
determined by random walk model and there is unknown 
linear equation which correctly describes of the transfer 
phenomenon of diffusion type. 

In this paper the simplest variant of the random walk 
model for ID system is considered. In this model the 
hoppings of a particle (or quasiparticles) are realized with 
equal probabilities during the interval of time At for dis¬ 
tances +Ax or —Ax . Then for the discrete variant of 
the problem formulated one can obtain the solution for 
the probability P of a particle location in coordinate x m 
in the moment t n in the form— 


P(%m 5 t n ) — 


(tn/At)! 


tn/At+x m /Ax , 2(i„/Ai) 


( 7 ) 

It is known that expression ( 2 ) presents one of the partial 
cases of solution of the random walk problem describing 
the probability density distribution of a random value. 
Really, if divide ( 7 ) on 2 Ax and put 


( A X ) 

- > D = const > 0 (8) 

at Ace —> 0, At — > 0, then one can obtain expression for 
the probability density function describing a probabil¬ 
ity of location of diffusive particle in the vicinity of the 
poinl— x in the form ( 2 ). In this case, simple calculations 
show that the velocity V of perturbations propagation 
accepts the infinite value (V = lim(Ax/At) -> oo at 
Ax —> 0, At —► 0 ) but the solution itself has automodel 
character with respect to variable x 2 /t. We should note 
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that ir£ for simplification of expression (7) containing 
factorials with the usage of the Stirling’s decomposition 
formula 

*! ” *' 'XPMVSFI {1 + + jdjj + 0(s -3 ) j(9) 

only the first term in the braces is kept. It allows to in¬ 
troduce a criterion for correct evaluation of this approach 
in the form of the following inequalities 

t/At » 1, (10a) 



x m « tnAx/At. (10b) 

The fundamental solution for telegraph type equa¬ 
tion (4) represents another limiting case of the random 
walk problem^. In this case lim(Ax/Af) = V < oo 
at Ax — > 0,Af —> 0. The selection of fundamental so¬ 
lutions as a ground for analysis of transfer processes is 
justified by the fact that for linear equations of mathe¬ 
matical physics the knowledge of a fundamental solution 
determines specific solutions of these equations as well, 
if the corresponding initial and boundary conditions are 
known. It is related to the fact that the correspond¬ 
ing Green functions can be expressed in the form of lin¬ 
ear combinations of fundamental solutions including their 
derivatives and integrals. Then at transition to continu¬ 
ous variables it is convenient to measure a current length 
in the unit of elementary length Ax and the current 
time in the units of At. We put also V = Ax/At = 1, 
D = 1/2, r = V 2 /D = At/2 = 1/2. 

Then the fundamental solution of equation 0 is writ¬ 
ten as 

f ’ G(l -‘ ) = 7 is exp( “^ ) (11) 

and fundamental solution for (4) as 

Pte(x, t ) = exp(-t) J 0 (Vf 2 - x 2 ^ 0 (f - |ar|) (12) 

and the probability density distribution function of a par¬ 
ticle located in the vicinity of the coordinate x and in the 
moment time t for model (5) and considered in the given 
paper as the fundamental solution of unknown transfer 
equations is written as 


FIG. 1. Comparison of the values and for the fixed values of 
time t=30 and t=100. 
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FIG. 2. The differences of solutions (11), (12) and (14) for 
the values t=30 and t=100. 


Here and below the symbol G (Gauss) indicates on equa¬ 
tion Eq. o, the symbol TE (Telegrapher equation) on 
equation (5), and symbol RW (Random Walk) on initial 
equation (7). Then we will use the following hypothesis: 
the true solution of evolution equation of the diffusion 
type presents the random walk solution in the form of 
fundamental solution (14) but presentations (2, 11) and 
(6, 12) present approximate approaches. Solution (14) 
corresponds to the unknown linear evolution equation. 


III. RESULTS 


PRw{x,t ) = 


(*)! 


2 ! (M'] ! 2* 


(13) 


Or, being presented in the form of the Euler’s gamma- 
functions it can be presented as 


PRw{x,t ) = 


T(t+ 1) 


r + 1 ) r 2 *+! 


(14) 


Let us consider the evolution of equations (11, 12 and 
14). On Figurel we demonstrate the values of the func¬ 
tion pnw end its gradient V puw for the fixed times t= 
30 and 100. On Figure2 we show the differences of the 
function p^w (14) from the functions pc (11) and pte 
(12) calculated at the same values of time. 

With increasing of the time the functions considered 
not changing their forms are decreased in their values and 
are stretched in the space. The differences between the 
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functions pc, Pte , Prw become small and are decreased 
with the time growing. 

Taking into account the influence of the second decom¬ 
position term figuring in the Stirling’s decomposition ex¬ 
pression (11) for factorials and the asymptotic behavior 
of the Bessel function at (z » 1) 


Io(z) 


exp(z) 

■\J2 tt z 


(l+ j- z + ■..), 


(15) 


it is easy to calculate the values of the functions con¬ 
sidered at the point x = 0 in different values of time at 
t » 1 . These values are equaled approximately 


PRw(0,t) S= 

~ V2^t {1 41 + - )> 

(16a) 

Pte{ 0, t ) R 

' v / 2^t (1+ 8t + - ) ’ 

(16b) 

1 

II 

o' 

o 

(16c) 


From these relationships it follows that the differences 
of the functions (16) are decreased with time as t ~ 3 / 2 , 
their relative differences as f _1 and their absolute values 
as t -1 / 2 . The corresponding calculations show that for 
other values of the current coordinate x the differences of 
the functions (11, 12 and 14) values are decreased with 
time as t -3 / 2 . Incidentally, the spatial stretching of the 
functions themselves at t » 1 is proportional to t 1 / 2 . 

Let us consider the gradient concentrations dynamics 
corresponding to solutions (11), (12) and (14). For the 
function p^.w in notations used the corresponding gradi¬ 
ent will look as 


V prw 


1 F(t — l)(H l ((f — x)/2) — T((f + x)/2)) 

2 2 t+1 r((f — x)/2)T((t + x)/2) 


where ip(x) = ^ln(r(a;)) , and for the functions pc, 
Pte , as 


Vp G 


x 1 
t y/2nt 



(18) 


and 


V pte = —xexp(—t)h (\/t 2 — x 2 ^ / (V t 2 - x 2 ) . 

(19) 

For all functions considered the gradient concentrations 
distributions at the fixed moments of time coincide with 
the functions that tend to zero at x —> oo and x —> 0 and 
having one maximum. For equation (18) at the fixed 
value of time the maximum value should be observed 


at x m = a/2 Dt and the value of function (18) in this 
point is Vp G (a; m ) = l/(\/87r ~eDt) , the velocity of the 
movement of the maximum point V G m and its halfwidth 
will be equaled to V G m (x m ) = \JD/2t ~ t 1 ^ 2 and 
A® ~ 1.6y/2Dt ~ t 1 / 2 , correspondingly. The calcula¬ 
tions show that the differences between maximum loca¬ 
tions points of functions pc , Pte Prw and their relative 
differences are decayed as f -1 and the relative changes 
of the maximum values of the corresponding gradients 
V pg(x, t), Vpte(x, t) and \7prw(x, t) for large values of 
t are equaled 

X 7 n m 

1 - f_RW_ ~ 0.89961 f _1 , (20a) 

^Pte 

X 7 n m 

1 - ~ 0.54298 r 1 . (20b) 

' G 

Let us consider now the relationship between flux and 
concentration gradient. For the Gaussian (2) the follow¬ 
ing relationship between flux and concentration gradient 
is valid 


i^ ex p(-f^) _ d 7 exp (~^) ,. 

2 Dt y/A^Dt dt J yfAnDt 

X* 

For the functions pte(x, 1) and PRw{x,t) the fluxes in 
the points x* (x * <Vt) are calculated as 

vt 

d f 

J T E(x*,t) = — pTE(x,t)dx, (22a) 


Vt 

d f 

jRw { x *, t ) — — / p RW (x,t)dx. (22b) 

X* 

The calculations lead to the following dependencies of 
the square differences of the functions pg{x, t),pTE{x, t), 
PRw{x,t ), their gradients and the corresponding fluxes 
at V = 1: 


Irw-te = 


\ 


j (Prw — PTE) 2 dx ~ t 1:272 


(23a) 


Irw-g = 


\ 


J (PRW ~ P G f dx ~ t 1255 


(23b) 
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Ite-g = 


\ 


X — L 

J ( PTE - Pg ) 2 dx ~ t 


-1.269 


(23c) 


IVRW-VTE = 


\ 


J (V prw — VpTE) 2 dx ~ t 1626 


(23d) 


IVRW-VG = 


\ 


J (V Prw - Vpcfdx ~ r 17335 


IVTE-VG = 


N 


(23e) 

X—t 

J (ypTE-Vpcfdx-t- 1 ™, (23f) 


LJrw—Jte 


\ 


x—z 

J (Jrw ~ Jte) 2 Jx ~ r 


1.022 


(23g) 


'-Jrw — Jg 


\ 


X—Z 

I (Jrw ~ Jcfdx ~ r 1 032 , (23h) 


1 -Jte — Jg 


\ 


x—z 

J (Jte — Jc) 2 dx ~ i^ 1 029 . (23i) 


One can notice that deviations of the values considered 
from each other are decreased with time faster than i _1 . 


IV. DISCUSSION 

So, for the temporal intervals t > 10 3 At calculating 
errors in replacement of the ’’accurate” solution by the 
functions , at calculations of the desired fluxes will be 
less then 10~ 3 and for the times t > 10 6 Af will be less 
then 10 -6 . We should remark also that absolute values of 
the fluxes and their concentration functions are relatively 
large only for the times x < 10 _1 t. For the times t = 10 2 
the flux value is becoming less then 10~ 3 at x = 20, for 
t = 10 3 is less then 10 -3 at x = 200, for t = 10 4 is less 
then 10 -3 at x = 2000. 

Let us consider the changings of deviations from 
Urw + xdjjiw/dt) for the function pnw(x, t). These de¬ 
viations are defined as 



FIG. 3. The values of the corrected function F at various 
values of the current coordinate x. 


L'^PRW — JRW — 'rdjRw/dt 


\ 


X — Z 

J (V prw — (Jrw + TdJRw/dt)~dx (24) 


It was proved that function (24) is decreased quickly with 
growth of the time as I v p RW -j RW -TdjR W /dt ~ t~ 2 - 
We introduce for equation (4) the corrected function 
F(x,t) by means of relationship 


Jrw + xdJuw j dt = —DF(x , t)VCuw(x, t) (25) 

Here the index RW stresses the origin of this relationship 
with random walk problem. Really, relationship (25) can 
present in the form of (4) for the random walk problem 
for ID case with replacement of the coefficient D for the 
parameter Drw = DF(x, t) . The analysis of (25) shows 
that with the increasing of the observation time F —> 1 , 
but for the various values of the coordinate x the times 
of approaching of the function F to the unit value lare 
different as it is shown on Fig.3. 

We want to note here that simple and good enough 
analytical approximation for the function F at the values 
x » 1 and t » 1 serves the following expression 


F(x,t ) ss (1 — exp(— 2 Vt/ |:r|))0(VT — |x|) (26) 

In this expression the function F is presented as the ap¬ 
proximation to the Heaviside function 0 for each fixed 
value of the coordinate x and depends only on the char¬ 
acteristic velocity V associated with random walk step. 
From this observation it is easy to conclude that for each 
moment of time t e if we define deviation of the function 
F from the unit value as e = 1 - F the value of the 
corresponding coordinate corresponding to the interval 
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t >> 1 can be written as x e ss (1/21?) (lne 1 ) t e . If 

we put, for example, e = 10 -3 , V = 1 then x e = 0.29 t e . 
It signifies that the relative velocity of approaching of F 
to the unit value (F —> 1) is considerably higher in com¬ 
parison with the propagation velocity of the fixed density 
distribution value proportional to t 1 / 2 . 

Then equation (25) can be rewritten in the form 


J + rdJjdt = —D( 1 — exp(— 2 Vt/ |x|)) 

VC{x,t) <d(Vt-\x\) (27) 

where we use the notations introduced above V = 
Ax/At, D = (Ax) 2 /(2At), t = At/2 , i.e. they are 
expressed by means of microparameters associated with 
random walk problem. 

We should note also the influence of fluctuations. Re¬ 
ally, it is necessary to take into account that in derivation 
of equation m the locality principle is used, as well. It 
imposes certain limitations of the spatial scales where the 
diffusion type process takes place. For the equilibrium 
values of concentration associated with number of par¬ 
ticles and temperature one can use evaluations of their 
fluctuations in correspondence with formulas 


\/(( ajv > ! } i #*> i 
N ~ V7T T ~ v^/V 


(28) 


The following question can be posed, what is the size of 
a system should be to take into account the thermody¬ 
namic fluctuations of different macroscopic parameters 
as particles, temperature and etc.? It is obvious that 
we cannot show the accurate limits and it has a relative 
character. For example, for subsystem having 10 6 parti¬ 
cles the relative fluctuations give the value ICR 3 . But for 
small values of time and distances the thermodynamic 
fluctuations become very large and it is impossible to de¬ 
scribe propagation perturbations front using only deter¬ 
ministic models. In other words, for the times t < 10 3 Af 
and distances x < 10 3 Aa; the usage of mesoscopic ap¬ 
proaches becomes necessary. In this case we go out from 
the principle locality limits and continuous medium ap¬ 
proach. This statement coincides with the conclusions 
given above associated with comparison of the functions 
m and (5) and ’’true” function (14). 


V. CONCLUSION 

In this paper we consider local-equilibrium and non¬ 
equilibrium models of transfer phenomena of the diffu¬ 


sion type. These phenomena are considered in the frame¬ 
work of the supposition that true solutions correspond to 
the solutions of the random walk problem (14) while the 
models (1-3) and (4-5) are considered as approximations 
to the ’’true” solution. The comparisons given above and 
analysis of the differences in fundamental solutions for 
equations (jT|), (5) and continuous representation of the 
random walk function as the density distribution func¬ 
tion for random value (14) allow to make the following 
conclusions: (a) the differences between the values of the 
functions considered are decreased with time as f -3 / 2 ; 
(b) the propagation velocity of the fixed values of the 
functions considered is proportional to t 1 / 2 ; (c) on rela¬ 
tively large values of observation times and distances the 
fundamental solutions becomes very close to each other. 

Taking into account these conclusions one can formu¬ 
late the criterion of applicability of the simplest equations 
as the telegraph equation (5) or heat transfer ([!]). Their 
criterion of applicability can be expressed in the form 
of inequalities: t/At >> 1, x « tAx/At and t » r 
where Ax, At, r are defined as microparameters. We 
should note that for the most problems of macrophysics 
that conditions are satisfied for relatively small values of 
the observation time t ~ 10 2 A< . For more accurate cal¬ 
culations of the desired relationships between fluxes and 
concentration gradients in the framework of random walk 
problem one can use relationships (25-27). In these rela¬ 
tionships the additional parameters figuring in diffusion 
equations are expressed in the form of microparameters 
related to the random walk problem. 
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